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We discuss extensions of a recently introduced model of multi-critical CDT to higher 



■^ ■ multi-critical points. As in the case of pure CDT the continuum limit can be taken on 

the level of the action and the resulting continuum surface model is again described by 
lO ■ a matrix model. The resolvent, a simple observable of the quantum geometry which is 

^f) \ accessible from the matrix model is calculated for arbitrary multi-critical points. We 

O ■ go beyond the matrix model by determining the propagator using the peeling procedure 

which is used to extract the effective quantum Hamiltonian and the fractal dimension in 
agreement with earlier results by Ambj0rn et al. With this at hand a string field theory 
formalism for multi-critical CDT is introduced and it is shown that the Dyson-Schwinger 
equations match the loop equations of the matrix model. We conclude by commenting on 
how to formally obtain the sum over topologies and a relation to stochastic quantisation. 
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1 Introduction 

A long standing problem in the causal dynamical triangulation (CDT) approach to causal 
quantum gravity (see [2] for reviews on recent progress in 2D and 4D) has been to obtain 
an analytical understanding of the behaviour of models that include matter. Shortly 
after the original proposal of CDT in [2], a number of models appeared generalising 
CDT which exhibited scaling behaviour different from pure CDT [3], thereby potentially 
corresponding to some form of matter coupled to CDT. One short-coming of this early 
work was that the models did not have a clear flat space analogue and therefore it was 
difficult to make a comparison to flat space theories or even to matter models coupled to 
dynamical triangulations (DT), i.e. two-dimensional Euclidean quantum gravity. For DT 
the known analytical results are for the case when the matter sector corresponds to one 
of the {p,q) minimal conformal field theories (CFTs). It is these CFT models that are 
of greatest interest, since they may be defined both on non-dynamical lattices and DT. 
However, until recently models of such CFTs coupled to CDT have resisted an analytical 
treatment. 

Recent progress in [U [5] has changed this situation. In |1], the matrix model formu- 
lation of CDT proposed in [6l [7] was generalised to the case of a quartic potential, and it 
was shown to possess a multi-critical point whose scaling exponents put it in a universality 
class different from pure CDT. In DT such a matrix model corresponds to hard dimers on 
a random lattice and therefore to a (2, 5) minimal model in the continuum limit. A priori 
the interpretation of the multi-critical point in the CDT model was less clear; indeed if 
one adds polygons of arbitrary type to a CDT in a transfer matrix formulation one does 
not obtain any new critical points. For dimers this can be understood as being due to 
the fact that each square corresponds to a dimer orientated in two different ways, only 
one of which respects the casual structure. However, it was argued in [4j that the new 
multi-critical point indeed corresponds to dimers coupled to CDT. This argument was 
confirmed by simultaneous work [5] in which a hard dimer model on a CDT was analysed 
directly via a combinatorial argument based on tree bijections. This givs us confidence 
that the multi-critical points of the CDT matrix model correspond to the (2, 2p+l) CFTs 
coupled to CDT, as is the case in DT. 

Given these developments, some obvious but nonetheless interesting and important 
quantities characterising the quantum geometry of multi-critical CDT remain to be cal- 
culated. In this paper we discuss generalisations of the work in [1] ; a CDT matrix model 
with an arbitrary potential possess a hierarchy of multi-critical points, the number of 
which is only limited by the order of the potential. We note that similarly to the case 
without any matter coupling |6], the CDT scaling limit may be taken at the level of the 
action, resulting in a matrix model formulation of the continuum surface model. At this 
point we would like to stress that this observation was already made in [7] . However, in [7] 
the resulting expressions were wrongly interpreted as belonging to the same universality 
class. To clarify this, we redo the analysis and explicitly determine the loop equations 
and resolvents for an arbitrary multi-critical point. 

Having this relation at hand allows for a whole host of previous results [H [9l [TOl [11] 



to be easily generalised which is the main content of this paper. In particular, as a gen- 
eralisation of [8] we determine the propagator using the peeling procedure in a manner 
similar to that used in [12] for DT. From this we extract the effective quantum Hamilto- 
nian and the two-point correlation function. Using the latter we compute the Hausdorff 
dimension du- We find in agreement with the conjecture in |1] that d^ = m/{m — 1) for 
the mth multi-critical point. Besides the grand-canonical Hasudorff dimension we also 
calculate the canonical Hausdorff dimension which agrees with the former. Having derived 
the effective quantum Hamiltonian, we set up a string field theory (SFT) formalism for 
multi-critical CDT, analogous to the case of CDT without matter coupling [9j and DT 
[T3] . It is shown that the Dyson-Schwinger equations (DSE) of the SFT reproduce the 
loop equation of the matrix model, which can be used to determine the loop amplitudes 
order by order in the topological expansion. We conclude by discussing possible exten- 
sions related to the sum over topologies and stochastic quantisation in a similar vein to 
what has been done fTOl HI] for CDT without matter. 



2 Multi-critical Points of the CDT Matrix Model 

Consider the one hermitian matrix model, 

Z = /"[#]e-f ^^(<^), (1) 

where V{(f)) = —gcj) + |0^ — |0^. It possesses the loop equation, 

mA^? = V\z)W^{z) + W^{z), (2) 

where W^</,(-2) is the resolvent for and W^{z) is a polynomial of order one. Note that since 
the resolvent corresponds to the disc amplitude with boundary cosmological constant z, 
it can be expressed as W^{z) = YlkLo z~^~^wi where wi is the disc amplitude for a disc 
with a boundary of fixed length /. 

It was first shown in [7j that ([2]) can be interpreted as describing the combinatorics 
of a DT in which spatial topology change i.e. the birth of a baby universe, is assigned a 
coupling constant (3. Explicitly this is because (|2]) may be rewritten for large enough / as, 

1-2 
Wi = gWi+i + gWi^i + /3 ^ Wi^m~-2Wm (3) 

m=0 

which can be interpreted graphically as shown in Figure [H Note that the linear term 
in the potential can be interpreted as the second term on the right hand side of Figure 
[H [7]; it adds a triangle (factor of g) and increases the boundary by one. Interpreting 
the loop equation combinatorially we see that if we do not include the linear term the 
boundary cannot increase once /3 is taken to zero. It is then clear that the term with 
coefficient /3 is the only one to allow spatial topology change. Hence, by taking (3 to 
zero the contribution from geometries in which spatial topology change occurs can be 





Figure 1: Graphical illustration of the loop equation of normal CDT. 



suppressed. Furthermore, it was shown that there exists a non-trivial scaling limit about 
the point /3 = corresponding to continuous geometries in which topology change is also 
assigned a scaled coupling constant Qs- It has long been known [2] that the key difference 
between DT and CDT is the absence of topology change in the latter and therefore this 
new scaling limit was interpreted as corresponding to a "generalised CDT" in which a 
controlled amount of topology change is permitted. This interpretation is justified given 
that known CDT results are recovered upon setting Qs = 0. 

A particularly interesting feature of the amplitudes obtained for generalised CDT 
is that they can be written in a form that matches the unsealed functional form of the 
amplitude, thereby suggesting that the continuum quantities can be obtained directly from 
a matrix model. That this is indeed the case was shown in [7] by demonstrating that the 
scaling limit associated with CDT could be taken at the level of the action. In particular 
by substituting into ([T]) the scaling forms, 

= 0* Inxn + a^, g = g* - a^K and /3 = a^gs, (4) 

where lAfxiv is the unity matrix and 0*, and g* are defined by requiring V'{g*,(l)*) = 
V"{g*,(l)*) = 0, we obtain, 

"[rf^Je-^^'-t^^^-i^l, (5) 

where Z is the continuum version of Z, obtained via a multiplicative renormalisation. The 
partition function in (^ defined the continuum theory of generalised CDT. Amazingly, 
since the size of the matrix A^ in ([S]) only plays the role of a coupling constant in the 
perturbative expansion, one can obtain a summation over the entire perturbative string 
expansion simply by setting A^ = 1 and evaluating the above integral [10]. 



In the remainder of this section, we now repeat the analysis of [7] in detail to generalise 
this above fact to all multi-critical points, including the first multi-critical point found 
in in [5], where in contrast to [7] we interpret the higher multi-critical points as new 
universality classes. One first extends the matrix model so that the potential has the 
form V^(0) = —g(j) + |0^ — g^^^^ ^'^^ ■ '^'^^^ leads to loop equations of the form (|2]), 
where W'^ is now of order ?7i — 1, and the corresponding length-space loop equation, 

m+l 1-2 

Wi = gWi^i +g^ tkWl+k-2 + /3 X] ^l-m-2Wm- (6) 

k=3 m=0 



Note that we let ts = 1, since in the combinatorial interpretation of the loop equations, 
shown in Figure [1] and generalised above, this corresponds to each triangle having a 
fugacity of g. We therefore have m — 1 free parameters in the potential. 

Although one could solve the loop equation and then look for a scaling limit of the form 
found for the pure model, it is instead quicker to realise we can again take the scaling 
limit about /3 = at the level of the action. To this end, consider the critical point 
({tfc},5f*,0*), it is defined by the condition that the derivatives V^'^\{tl],g* ^cjf) = for 
1 < n < m. This constitutes m equations for the m — 1 parameters in the potential 
together with the quantity (j)*; barring any degeneracy we expect to be able to solve for 
the critical point. We then substitute the scaling forms, 

= 0* l^x^ + a$, g = g* - a^A and /3 = a^'^^gs, (7) 

into the action to obtain, 



where we have introduced the notation V* for V evaluated at the critical point. Again 
we performed a wave function renormalisation by absorbing multiplicative factors in Z. 
In general, to obtain a non-trivial scaling limit the critical point introduced above should 
correspond to a point in the parameter space of Z at which Z is non-analytic. For example, 
in the case of DT the scaling limit corresponds to zooming in on the region near a branch 
point of Z. Here, by taking the scaling limit at the level of the action we have made 
this point highly implicit, indeed a priori it is not clear the scaling limit defined here will 
produce a non-trivial result. However, it is clear the integral in (|8]) resulting from the 
scaling limit is non-trivial. Furthermore, since we do not scale A^ in this limit we can see 
that the scaling limit actually commutes with the matrix integral in ([T]). To be explicit, 
([T]) may be written in terms of the eigenvalues Aj of 0, thus, 

N 

\{dX.\{{\-\,fe-Ty^'^\ (9) 

The scaling ansatz for Aj is, from ([7]), Aj = 0* + aAj. Subbing this into ([9]) together 
with ([7]) we see that the measure is only changed by an overall factor of a. We therefore 
conclude that the small a expansion of ([9]) can be found either by computing the small 
a expansion of the integrand and then performing the A^ integrations of the eigenvalues, 
as we do in this article, or by first computing the integrals over the eigenvalues and then 
taking the small a expansion of the result as done in [7], |4]. The advantage of doing 
the latter is that we see explicitly that this scaling limit corresponds to zooming in on a 
branchpoint of the partition function (for a disc) while simultaneously taking the branch 
cut length to zero. Here we are content with commuting the a — )■ limit through the 
integral and viewing the resulting non-trivial integrand as evidence that the critical point 
is non-trivial. Furthermore we see that any scaled amplitudes may be computed from the 
matrix model (jH]) without any scaling. 



Let us for simplicity choose the coefficients in the action of ([H]) such that the potential 
becomes 

V(X) = AX + ^^(-Xr+\ (10) 

m + 1 

One can now easily derive the loop equation for the resolvent W{X) of the matrix model 
with the above potential, 

g,W{Xy = [A - {-X)'^]W{X) + W^{X), (11) 

where W'^lX) = —{—X)"^^^+Yy^ aiX^ is a polynomial of order ?7i—l and a; = ai{A,gs). 
The solution of the loop equation (ITT!) is formally given by 

^^^) ^ -((-X)- - A) + v/((-X)"- - Ay + Ag,WQ{X) ^^^^ 

■^9s 



where W^{X) is still undetermined. However, one can fix W^{X) by noting that (TT2l) 
has to obey the one-cut solution 

^(^) ^ -((-X)- - A) + (x-^ + Er=o' fe^x V(X - cO(X - c) ^ ^^g^ 



where the 6;, ci and C2 are functions of A and Qs- Matching (1121) and ( 1131) yields a set of 2m 
equations for the 2m unknowns a;, 6^, Ci and C2 which completely fix the solution. Solving 
the set of equations can be tedious for general m and (7^ > and we will not peruse it 
here. However, we note that in the case of most physical interest of pure multi-critical 
CDT, i.e. when Qs = 0, the resolvent simplifies to 

W{X) = , ^J}^^ , . (14) 

Demanding that the inverse Laplace transform W{L) — ;■ as L — )■ 00 and furthermore 
that there exist no oscillatory terms in the large L limit jj we find, 

m— 1 

W'^iX) = Y^i-iy+'A^'^X^. (15) 

This yields for the resolvent 

WiX) = —-^——. (16) 

In agreement with the corresponding result for m = 3 as obtained in |3]. Furthermore, it 
can also be shown that the combinatorial approach derived for (restricted) dimers coupled 
to CDT [5] yields the same result with m = 3. 



^This requirement is equivalent to the one-cut assumption. 



3 Finite Time Propagator and Fractal Dimension 

Although the matrix model approach is well adapted to analytically studying different 
forms of matter coupled to CDT it is not so suited to providing information about the 
geometry of the spacetime. An important observable encoding more refined properties of 
the geometry is the finite time propagator G{x,y;t), which is unavailable in the matrix 
model. The quantity G{x,y;t) corresponds to a summation over manifolds of cylindrical 
topology in which the entrance and exit loops are separated by a fixed geodesic distance 
t and with the boundaries carrying the boundary cosmological constants x and y respec- 
tively. The finite time propagator is also a necessary ingredient in the construction of the 
string field theory which is the topic of the next section. In this section we will compute 
G{x, y; t) and discuss how a naive calculation of the fractal dimension dn agrees with the 
value conjectured in [4j from considerations of multi-critical branched polymers. 

To begin our calculation we first introduce G{1, /'; t), which is the length-space version 
of the propagator, in which the boundaries are of fixed length / and /'. We then have 
G{x,y;t) = Ylu=oYll/'=o^~''~^y~'' ~^^i^^^''^'^)- The finite time propagator may then be 
computed, for the case of a general potential in the matrix model, by viewing ([6]) as a 
time dependent process, as first suggested in [12]. We first introduce a new boundary into 
the spacetime, leading to a modification of ([6]) of the form, 

m+l 1-2 

G{1, n = gG{l - 1, /') + ^ 5^ hG{l + A; - 2, r) + 2/3 5^ wi^m-2G{m, I'), (17) 

fc=3 m=0 

where G{1, 1') is the propagator with no constraint on the distance between the initial 
and final boundary and the factor of two in front of the topology changing term is due 
to the ability to add the exit loop to either of the two baby-universes. Now we envision 
the addition of polygons to this amplitude as a time dependent process, this is known 
as the peeling procedure. The LHS of the above equation has one additional polygon on 
its boundary compared to the RHS. If we were to add an entire space-like layer of new 
polygons to the amplitude, the distance between the initial and final boundary would be 
increased by one. Hence for large enough / we can approximate the addition of a single 
polygon by a change in t of 1//, leading to, 

\dtG{l,l';t) = gG{l-l,l';t)-G{l,l';t) 



I 



m+l 1-2 



+g Y^ tkG{l + k-2, /'; t) + 2(3j2 Wi-m-2G{m, I'; t). (18) 

fc=3 m=0 

By summing over / and I', this can then be written as, 

dtG{x, y; t) = d, {{V'{x) - 2(3W^{x))G{x, y; t)) . (19) 

Substituting in the scaling relations ([7]) one find that to obtain a non-trivial limit we must 
have t scaling as t = a^~^T. After rescaling A and X we obtain, 

dTG{X, Y- T) = ~dx [((-X)" - A + 2gM{X)) G(X, F; T)] . (20) 



Note that G{X, Y; T) and the corresponding expression in length space are related by a 
Laplace transformation 



/•oo noo 

Jo Jo 



(21) 



To solve the differential equation f l2U]) one has to set the initial conditions. In particular, 
one imposes that G{Li, L2; 0) has only support on Li = L2. This implies that G{Li, L2; 0) 
is given by a sum over 5^^\Li — L2), n > 0, where S^"'\-) is the n-th derivative of the 
delta function. In the Laplace transform this is equivalent to 

G{X,Y;0) = p^, (22) 

where P(X) is a polynomial in X. We will now show that consistency with the results 
of the previous section requires that P{X) = dxW^{X) with W'^lX) given in (llSp . In 
particular, we show that (1201) together with (I22p yields the correct disc function. 

Let us first note that the disc function satisfies W{X) = J G{X, L2 = 0; T) dT with 
2gs — !■ Qs, i.e. the disc function corresponds to a propagator of arbitrary time with the final 
boundary shrunken to zeroo Integrating fl20|) with respect to T and shrinking L2 — )■ 0, 
or equivalently F — t- 00, and removing the symmetry factor 2 one obtains using fl22l) that 

- dxW'^iX) = -dx [ii-Xr - A + 9sWiX)) WiX)] (23) 

which after integrating with respect to X precisely yields fITT]) . 

One observes that (1201) yields the following Hamiltonian for pure multi-critical CDT, 
i.e. with Qs = 0, 

- drGoiX, Y- T) = Ho{X)Go{X, Y; T), Ho{X) = ^^((-X)"' - A), (24) 

Using the fact that G{X, Y; T) and G{Li, L2] T) are related by an inverse Laplace trans- 
formation with respect to X and Y one obtains the Hamiltonian in length space 

Ho{L) = -Ld'^ + AL. (25) 

Note that for m even the Hamiltonian is Hermitian with respect to the measure dfi{L) = 
L-^dL. 

We now turn our attention to the analysis of the fractal structure of the space-time 
ensemble as characterised by its fractal or Hausdorff dimension. Given that for any value 
of Qs, the production of baby-universes is under control, we expect the fractal structure 
of such space-times to be equivalent to the f^^ = case. We therefore focus our attention 
on the case of pure multi-critical CDT, by setting gg = 0. 



^Note that removing the symmetry factor 2gs — )■ gs is the same as in pure CDT ^9 and similarly in 
the DT string field theory [13] . It is naturally, keeping in mind that the propagator is not the two- loop 
function, as is for example explained in |^. 



Using the method of characteristics to solve the above equation gives, 

Go(X, y; T) ^ dxW<^iXiT,X))i-XiT,X)r-A 
^ ^ X{T,X) + Y (-X)™-A ^ ^ 



with 

dX{T,X) 
dT 



-((-X(T,X)r-A), (27) 



and the condition X(0) = X. One can now shrink the initial and final boundary to zero 
by taking X and Y to infinity, resulting in the amplitude Gq{T) for a sphere with two 
punctures separated by a distance T, 

G'o(T) = axVr^(X(T)) [(-X(T))™ - AJ , (28) 

where we set X{T) = X{T,X = oo). One now has 

From this expression one already sees that one has the following scaling 

Go(T) = A2^^Fi(TA^^), (30) 

in terms of a yet underdetermined function Fi. This scaling shows that area scales as 
A ~ j'm/im-i) ^^^ thus that the Hausdorff dimension is 

dn = 31 

m — 1 

in agreement with the value of [3] suggested from the comparison with multi-critical 
branched polymers |14j , as well with a recently introduced model of multi-critical tensor 
models [15] . This dimension can be seen as a "grand-canonical" definition of the Hausdorff 
dimension [16]. One can also define the "canonical" Hausdorff dimension dh- To do 
so consider the volume of a spherical shell {L{T))y at fixed volume V ^ T'^^. The 
"canonical" Hausdorff dimension is then defined through {L{T))y ~ T'^^~^ . In particular, 
one expects the following scaling [THl E] 

{L{T))y = V^-^I^"F2{TIV^'''^), (32) 

where F2{x) ~ a;^h-i g^g j; _>. q g^j^fj ^2(0;) ~ exp(— a;'^^/^"'^^^'') as x — t- 00. For pure 
CDT one has for example dn = dh = 2 and for DT dn = dh = 4. However, for multi- 
critical branched polymers it is known that the both dimensions differ and one has du = 
m/{m — 1), while dh = 2 [T^ [T7] . We have already seen that for multi-critical CDT we 
have dn = m/{m — 1) in agreement with the multi-critical branched polymers. In this 
comparison it is therefore interesting to see whether for multi-critical CDT one has the 



same behaviour. We can obtain the Laplace transform of {L{T))y from the propagator 
(!26|l . In particular, taking the first boundary to zero {X — )■ oo), one gets 



{Lin, 



G{Y;T)dY 



d 

G{Y;T) 



= X{T)-\ (33) 

Y=0 



For small T ^ l/A^/'^^ one has the following expansion from the characteristic equation 

(EHD, 



m + l 



(L(T))^ ~ T— - CmAT— + ... (34) 

However, from the Laplace transform of (1321) one expects that {L{T))j^ ~ j^dh-i ^qj- 
T <^ 1/A^/^^ and thus we conclude that dh = dn = ^;^- This means that the continuum 
model has the same fractal properties at all scales. This is probably what one would expect 
for a continuum surface model which are not simply trees as in the case of multi-critical 
branched polymers. 



4 String Field Theory 

In this section we introduce a string field theory (SFT) formalism for multi-critical CDT 
analogous to the SFT for generalised CDT [9], i.e. m = 2. A previous construction for 
DT can be found in [13]. 



Let us start by defining a vacuum |0) and creation and annihilation operators with 
the following properties 

[^(L),*t(L')] = M(L-L'), ^(L)|0) = (0|^^(L) = 0. (35) 

One then defines the "second-quantised" Hamiltonian 

— ¥{L)Ho{L)^iL), (36) 

where Hq{L) is the first-quantised Hamiltonian of pure mult i- critical CDT as given in 
(l25l) . From this one can rewrite the propagator G(Li, L2; T) for pure multi-critical CDT 

as 

Go(Li,L2;T) = (0|^(L2)e-^^»^t(^^)|o), (37) 

i.e. as the creation of a string of length Li, a propagation over time T with respect to 
Hq{L) and the annihilation of the string given that it has length L2. 

We can now set up the SFT by allowing strings to split, merge and be annihilated. 
This is implemented using the following Hamiltonian 

H = Ho-gsfdL, I dL2'^\L,)^^{L2)^{L, + L^) 

-agsjdL.j dL2^\U + L^)^ {L^)-^ {L{) - j^ piL)^iL), (38) 

10 
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L,+L 




(c) 



(d) 



Figure 2: The Feyninan rules of the string field theory: (a) represents the propagation of 
a single string with no spatial topology change, (b) and (c) correspond to the splitting 
and merging of strings respectively and (d) is a tadpole diagram representing the process 
of a string disappearing into the vacuum. 



as illustrated in Figure HI We introduced the coupling gs for each splitting and joining of 
the string as well as a factor of a controlling the joining process and thus the topology 
of the string world-sheet. In relation to the matrix model formulation it will be clear 
that a = 1/iV^, i.e. a = corresponds to the planar or large-A^ limit. Finally, p{L) is 
the tadpole term related to the annihilation of a string. One expects p{L) to only have 
support on L = and furthermore it is natural to have 

p(L) = G'(L,0;0), (39) 

where G{L, L2; T) is the propagator as derived in the previous section. 

Given the STF Hamiltonian fl38p . one can now use the field theoretical formalism (see 
[HI [13]) to derive n-loop amplitudes by introducing a source term in the Hamiltonian 
yielding 

(5"Z(J) 



lim ( 

r-i>oo 



|e-^^v]>t(Li)---vI/t(L„)|0) 



6J{Li)---6J{Ln) 
where the generating function Z{J) is given by 

Z{J)= lim(0|e-^^e^^^"(^)*'(^)|0). 



j=o 



T-!>oo 



(40) 



(41) 



The amplitudes limT-s.oo(0|e~"^''^ \l/^(Li) ■ ■ • \1'"''(L„)|0) still includes disconnected world- 
sheets. The amplitudes for the connected world-sheets are obtained from the generating 
function 



in a similar manner 

W{L,,...,L^) 



F{J)= log Z{J) 
«!-" lim (0|e-^^ ^!\Li) ■ ■ ■ ^^(L,)|0), 



(42) 



r-!>oo 



a 



l-n 



6"F{J) 



6J{Li)---6J{Lr, 



(43) 



j=o 
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Here W{Li, ..., L„) is the so-called n-loop amplitude. Note that we introduced the coupling 
a in the SFT Hamiltonian fl38|) precisely to keep track of the topology of the string world- 
sheet. However, a genus zero n-loop amplitude naturally comes with n — 1 mergers of 
the string. To compensate for this we introduced the factor of a;^~" in the definition of 
Ty(Li,...,L„). 

It was shown in [9] that commuting the \&(L)'s in fHT]) past the source term yields the 
following equation for F{J), 



6^F{J) 






(44) 



Taking derivatives with respect to 6/6 J{L) of this equation yields the Dyson- Schwinger 
equations (DSE) of the SFT. In particular, taking the first derivative with respect to J{L), 
setting J{L) = and Laplace transforming we obtain the DSE for the disc function, 

= Ho{X)W{X) - dxW^{X) + gsdx{aW{X,X) + W{X)W{X)). (45) 

Setting a = 0, i.e. in the planar limit, this equation precisely yields the loop equation 
(llip derived from the matrix model with potential (llOp in the planar limit. Furthermore, 
including the 1/N'^ contribution in the loop equation, fITT]) becomes 

gs{W{Xf + -^WiX, X)) = V'{X)W{X) + W'^{X). (46) 



Hence, differentiating the loop equation f l46|) with respect to X one obtains precise agree- 
ment with the DSE equation fj45l) if we identify a = l/N"^, as was also the case for the 
SFT of CDT without matter coupling [5]. Furthermore, one can convince oneself that 
the same equivalence between the DSE of the SFT on the one hand and the loop equa- 
tions of the matrix model on the other hand still holds true when considering higher loop 
amplitudes. 



5 Discussion 

In this paper we discussed how the multi-critical CDT scaling limit of the one matrix model 
found in [4] can be generalised to obtain a hierarchy of multi-critical points completely 
analogous to the multi-critical points of the DT scaling limit as obtained in [H [7]. In the 
DT case these mult i- critical points correspond to a continuum gravity model coupled to 
minimal CFT matter of (2, 2p -|- 1) type. Based on the analysis in |4l, i5j it is highly likely 
that a similar interpretation can be given to the higher order critical points. 

In going beyond the matrix model, we have obtained via a peeling produce the dif- 
ferential equation satisfied by the finite time propagator and hence the effective quantum 
Hamiltonian for the time evolution of the space- like universe. This allowed us to generalise 
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many of the results for generalised CDT without matter, already present in the literature. 
Firstly we have investigated the conjecture in |4j that the Hausdorff dimension (fractal 
dimension) of the multi-critical CDT is dn = m/{m — 1) both in the grand-canonical as 
well as canonical sense. Secondly, we have shown that a string field theory can be built 
using this Hamiltonian, which in principle allows for the calculation of any amplitude with 
a specified number of baby-universes and handles. The structure of this SFT is similar to 
the SFT constructed to describe the multi-critical non-critical string. However, a pleasing 
aspect of the construction here is that the tadpole term can be naturally identified with 
the finite time propagator between a universe of finite size and of zero size over an interval 
of zero time. This is in contrast with the DT case in which the tadpole term had to be 
selected by hand. 

An important quality of the multi-critical CDT scaling limit is that the limit may be 
taken at the level of the matrix integral. That this can be done is particularly interest- 
ing, as the Dyson-Schwinger equations can be identified with the loop equations of the 
"continuum" matrix model provided we identify the coupling a appearing in the SFT for 
the merging of two universes with l/N"^, where A^ is the size of the matrix in the scaled 
matrix model. This is an extremely important point as it means the Feynman diagrams 
appearing in the large N expansion of the continuum matrix model have a direct inter- 
pretation as processes in the SFT. The only role of A^ in the continuum matrix model is 
to modify the coupling associated to the merging of two universes, thereby providing an 
independent weight associated to adding a handle to the worldsheet. If we do not care to 
distinguish between the processes of merging and splitting of universes then we are free 
to set A^ = 1 in the the continuum matrix integral. This is particularly profound as it 
reduces the matrix integral to a simple integral of one variable. We therefore have that 
the entire non-perturbative expression for the continuum partition function is, 

Z{A)= /"dze-^(^^+^(-^)"^'\ (47) 



That this is indeed the case can likely be obtained by considering CDT in the framework 
of stochastic quantisation, which we expect to yield the generalised Hamiltonian, 

H{L) = -Ld'j: + KL-g,L\ (48) 

This would generalise the Hamiltonian, which includes topology change, found for the case 
of pure CDT using this approach in fTTj. The solutions to the Wheeler-de-Witt equation 
in that case do indeed correspond to the expression obtained by evaluating f l47|) for the 
case m = 2. This is something we hope to pursue in future work. 

Finally, it is interesting to note the connections between the work here and the results 
of [18]. In |T8] the continuum limit of the exact partition function for a stack of A^ FZZT 
branes in the {1,171) background was computed and was found to be given by, 

Z= /"[d$]e-^4^*-;;^*'"^^], (49) 

where $ is an A^ x A^ hermitian matrix and A is a diagonal matrix whose entries are 
the boundary cosmological constants of each brane. The continuum limit of mth order 
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multi-critical CDT is therefore mathematically equivalent to the theory of a stack of A^ 
FZZT branes in a {l,m) background with each brane carrying an identical boundary 
cosmological constant. This also relates to a similar observation made in [1] for the case 
of pure CDT. The physical origin of this picture is certainly something that should be 
better understood. 
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